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A metallic electron system near an orientational symmetry breaking Pomeranchuk 
instability is characterized by a "soft" Fermi surface with enhanced collective fluc- 
tuations. We analyze fluctuation effects in a two-dimensional electron system on a 
square lattice in the vicinity of a Pomeranchuk instability with d-wave symmetry, 
using a phenomenological model which includes interactions with a small momentum 
transfer only. We compute the dynamical density correlations with a d-wave form 
factor for small momenta and frequencies, the dynamical effective interaction due to 
fluctuation exchange, and the electron self-energy. At the quantum critical point the 
density correlations and the dynamical forward scattering interaction diverge with a 
dynamical exponent z = 3. The singular forward scattering leads to large self-energy 
corrections, which destroy Fermi liquid behavior over the whole Fermi surface except 
near the Brillouin zone diagonal. The decay rate of single-particle excitations, which 
is related to the width of the peaks in the spectral function, exceeds the excitation 
energy in the low-energy limit. The dispersion of maxima in the spectra flattens 
strongly near those portions of the Fermi surface which are remote from the zone 
diagonal. 

PACS: 71.10.Fd, 71.18.+y, 74.20.Mn 

I. INTRODUCTION 

Under suitable circumstances electron-electron interactions can generate a spontaneous 
breaking of the rotation symmetry of an itinerant electron system without breaking trans- 
lation invariance. From a Fermi liquid perspective such an instability is driven by forward 
scattering interactions and leads to a symmetry breaking deformation of the Fermi surface. 
Alluding to a stability limit for forward scattering derived long ago by Pomeranchuk,— it 
is therefore frequently referred to as "Pomeranchuk instability" . On a square lattice, a 
Pomeranchuk instability with d x 2_ y 2 symmetry, where the Fermi surface expands along 
the k x axis and shrinks along the k y axis (or vice versa), was first considered for the t-J, 2 
Hubbard, 3,4 and extended Hubbard model.— Fermi surface symmetry breaking in fully 
isotropic (not lattice) systems has also been analyzed.^ A Pomeranchuk instability leads 
to a state with the same symmetry reduction as the "nematic" electron liquid defined 
by Kivelson et al.— in their discussion of similarities between doped Mott insulators with 
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charge stripe correlations and liquid crystal phases. 

A Pomeranchuk instability usually has to compete with other instabilities, but can 
also coexist with other symmetry breaking order. For example, in the two dimensional 
Hubbard model with a sizable next-to-nearest neighbor hopping and an electron density 
near van Hove filling a superconducting state with a d-wave deformed Fermi surface is 
stabilized at least at weak coupling. 8 Superconducting nematic states have also been 
included in a general classification of possible symmetry breaking patterns by Vojta et 
al.— In this work we will however focus on symmetry breaking Fermi surface deformations 
in an otherwise normal state. Since a Pomeranchuk instability of electrons on a lattice 
breaks only a discrete symmetry, no Goldstone mode exists and symmetry-broken states 
are stable also at finite temperature in d > 2 dimensions. 

Electron systems in the vicinity of a Pomeranchuk instability have peculiar properties 
due to a "soft" Fermi surface, which can be easily deformed by anisotropic perturbations. 
In particular, it was shown that dynamical fluctuations of such a soft Fermi surface lead 
to a strongly enhanced decay rate for single-particle excitations and thus to non-Fermi 
liquid behavior.— Close to a d x 2_ y 2 Pomeranchuk instability the decay rate is maximal 
near the k x and k y axes and minimal near the diagonal of the Brillouin zone. The decay 
rate was computed in random phase approximation (RPA) for a phenomenological model 
describing electrons on a square lattice which interact exclusively via scattering processes 
with small momentum transfers.— Subsequently it was shown that the putative continuous 
Pomeranchuk transition in this model is usually preempted by a first order transition 
at low temperatures, that is the Fermi surface symmetry changes abruptly before the 
fluctuations become truely critical . 11112 However, for reasonable choices of hopping and 
interaction parameters the system is nevertheless characterized by a drastically softened 
Fermi surface on the symmetric side of the first order transition, and hence by strongly 
enhanced Fermi surface fluctuations.^ Furthermore, adding a uniform repulsion to the 
forward scattering interaction, the first order transition can be suppressed, and for a 
favorable but not unphysical choice of model parameters even a genuine quantum critical 
point can be realized.-^ 

In this work we present a detailed analysis of dynamical Fermi surface fluctuations in 
the vicinity of a Pomeranchuk instability with d a .2_ 2/ 2-wave symmetry in two dimensions 
and their effect on single electron excitations. We compute the dynamical density correla- 
tions with a d-wave form factor for small momenta and frequencies, the dynamical effective 
interaction due to fluctuation exchange, and the electron self-energy, from which the spec- 
tral function for single-particle excitations is obtained. In the quantum critical regime 
the density correlations and the dynamical forward scattering interaction diverge with a 
singularity familiar from other quantum phase transitions in itinerant electron systems 
with dynamical exponent z = 3. The singular forward scattering leads to large self-energy 
contributions E(k, u), which are proportional to d\, where d k is a form factor with d-wave 
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symmetry. At the quantum critical point, ImS(ki?, to) scales to zero as |c<j| 2//3 . Fermi liquid 
behavior is thus destroyed over the whole Fermi surface except near the Brillouin zone 
diagonal. We also compute the crossover from non- Fermi to Fermi liquid behavior in case 
of a large finite correlation length £ in the ground state. In the quantum critical regime 
at low finite temperatures the self-energy has a part due to quantum fluctuations, which 
obeys (u /T) -scaling, but also a classical part, which is proportional to T£(T) near the 
Fermi surface. The classical contribution violates (u;/T)-scaling and was overlooked in 
Ref.0 

The article is structured as follows. In Sec. II we define our phenomenological model 
and review its mean-field phase diagram. In Sec. Ill we analyze the dynamical effective 
interaction of the model in a regime where the nearby Pomeranchuk instability leads to 
strong Fermi surface fluctuations. The low energy behavior of the self-energy in the pres- 
ence of these fluctuations is then computed in Sec. IV, and results for the corresponding 
spectral function are discussed in Sec. V. We finally summarize and discuss the possible 
role of d-wave Fermi surface fluctuations in cuprate superconductors in Sec. VI. 



II. F-MODEL 

To analyze fluctuation effects in the vicinity of a Pomeranchuk instability we consider a 
phenomenological lattice model with an interaction which drives a Fermi surface symmetry 
breaking, but no other instability. The model Hamiltonian, which has been introduced 
already in Ref . llC| reads 

H = e k n ka + ^- /kk'(q)n k (q)n k /(-q) , (1) 

k.cr k,k',q 

where e k is a single-particle dispersion, n k (q) = J2 a c k- q /2 o- c k+q/2,cr, and V is the volume 
of the system. Since the Pomeranchuk instability is driven by interactions with vanish- 
ing momentum transfers, that is forward scattering, we choose a function /kk'fa) which 
contributes only for relatively small momenta q, and refer to the above model as the 
"f-modeP . Clearly this model is adequate only if the Pomeranchuk instability dominates 
over other instabilities and fluctuations in the system. Otherwise it would have to be 
supplemented by interactions with large momentum transfers. 

For hopping amplitudes t, t', and t" between nearest, next-nearest, and third-nearest 
neighbors on a square lattice, respectively, the dispersion relation is given by 

e k = —2 [t(cos k x + cos k y ) + 2t' cos k x cos k y + t"(cos 2k x + cos Iky)} . (2) 

For a simplified treatment, which however fully captures the crucial physics, we consider 
an interaction of the form 



/kk'(q) = «(q) + g(q) <4dk' 



(3) 
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with u(q) > and g(q) < 0, and a form factor d^ with d x 2_ y 2 symmetry, such as d^ = 
cosk x — cosfc y . The coupling functions u(q) and g(q) vanish if |q| exceeds a certain 
small momentum cutoff A. This ansatz mimics the effective interaction in the forward 
scattering channel as obtained from renormalization group calculations 3 - and perturbation 
theory^ for the two-dimensional Hubbard model near van Hove filling. The uniform term 
originates directly from the repulsion between electrons and suppresses the electronic 
compressibility of the system. The d-wave term drives the Pomeranchuk instability. 

The mean-field solution of the f-model has been analyzed for various choices of param- 
eters in a series of recent articles . 11 ' 12 ^ 3 In the plane spanned by the chemical potential 
/i and temperature T the symmetry-broken phase is formed below a dome-shaped tran- 
sition line T c (/x) with a maximal transition temperature near van Hove filling. The phase 
transition is usually first order near the edges of the transition line, that is where T c is 
relatively low, and always second order near its center M> The two tricritical points at the 
ends of the second order transition line can be shifted to lower temperatures by a sizable 
uniform repulsion u included in /kk' ^ For a favorable choice of hopping and interaction 
parameters, with a finite t" and u > 0, one of the first order edges is suppressed com- 
pletely such that a quantum critical point is realized. Although quantum critical points 
are usually prevented by first order transitions at low temperatures in the f-model, the 
Fermi surface is nevertheless already very soft near the transition, such that fluctuations 
can be expected to be important.— 



III. EFFECTIVE INTERACTION 



In this section we derive and analyze the dynamical effective interaction for the f-model, 
which is closely related to the dynamical d-wave density correlation function^ 

POO 

N d (q,u) = -i / dte wt (Mq,tW-q,0)]> , (4) 
Jo 

where the d-wave density fluctuation operator is defined as 

nd(q) = ^dk^k(q) . (5) 

k 

and nd(q, t) is the corresponding dynamical operator in the Heisenberg picture. We first 
analyze A<d(q, v) and the effective interaction within RPA, and then discuss higher order 
corrections. 

The RPA result for the d-wave density correlation function in the f-model is simply 

N d (q,u) = - . (6) 

1 - 0(q) ng(q, v) 
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with the bare d-wave polarization function 

n2(g, V ) = -(^L /(£p+q/2 " £ " /(W2 " f } 4 • (7) 

; y (2tt) 2 zy + z0+-(e p+q/2 -e p _ q/2 ) p ^ 

The infinitesimal imaginary part in the denominator specifies that we consider retarded 
functions. Note that the coupling w(q) does not enter here because mixed polarization 
functions with a constant and a d-wave vertex vanish for small q. 

The RPA effective interaction is defined by the series of chain diagrams sketched in 
Fig. 1, yielding 

1 - «(q) n°(q, i/) 1 - (/(q) Il3(q, i/) 

where n°(q, z/) is the conventional polarization function, defined as 11° (q, z/) but without 
the form factor d^. 

Near the Pomeranchuk instability the denominator 1— g(q) n°(q, v) in Eqs. (6) and (8) 
becomes very small for q — > and — > 0, if f vanishes faster than q, while 1— it(q) Il (q, v) 
remains of order one or even larger. The effective interaction is then dominated by the 
second term in Eq. (8) and can be written as 

T k k' (q, v) = 9 S d (q, v) d k dy (9) 
with g = g(0) and the "dynamical Stoner factor" 

S d(^ v) = 1 TTFFoT V • ( 10 ) 

i -p(q) n d(q ; v) 

The d-wave density correlation function iV^q, v) is also proportional to Sd(q, v). In case 
of a second order transition the static Stoner factor 

S d = lim lim 5 d (q, v) (11) 

q^O v— >0 

diverges on the transition line. Near the first order transition obtained typically for low 
temperatures in the mean-field solution of the f-model, Sd is still drastically enhanced, 
that is of order ten and more.— 

Within RPA, the dynamical d-wave density fluctuations near the Pomeranchuk insta- 
bility and the corresponding singularity of the effective interaction are obviously deter- 
mined by the asymptotic behavior of the d-wave polarization function 11° (q, v) for small 
q and u, which we describe in the following; for derivations, see Appendix A. Analogous 
formulae for the expansion of the conventional polarization function II (q,^) in three 
dimensions have been derived already long ago in the context of almost ferromagnetic 
metals by Moriya.— 
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At zero frequency 11° is a real function of q, which can be expanded as 

n°(q,0)=a(T) + c(T)|q| 2 + O(|q| 4 ) (12) 
for small q. The coefficient a(T) is always negative and can be written as 

a(T) = J defie-^N^e) (13) 

with a weighted density of states N ( p(e) = J 8(e — e p ) dp. A low temperature (Som- 
merfeld) expansion yields 

a(T) = -NM ~ y N^(jm) T 2 + G(T 4 ) , (14) 

provided that N d 2(fx) and its second derivative are finite. The coefficient c(T) can be 
expressed as 



c(T) 



J de f'{e-n) 



-^W£„ a (e) - -^N' d2Ae (e) 



(15) 



where N d 2 At (e) = j 5(e — e p ) d p Ae p with the Laplacian A = dp x + dp y , and N d 2 v 2(e) = 

I (I^F ~~ e p) dpV 2 with v p = |Ve p |. Primes denote derivatives with respect to e. 
The sign of c(T) depends on the dispersion and other model parameters. In the low 
temperature limit it will be sufficient to use 

c = c(0) = -^N' d 2 Ae (ri - 1 JV^) (16) 

For finite frequencies n°(q, v) is generally complex. For small q and v with v/\cj\ — > its 
imaginary part behaves as 

Imn°(q,z/)^-p(q,T)^, (17) 

where p(q, T) > depends only on the direction q of q. At low temperatures it is sufficient 
to use p(q) = p(q, 0), which can be expressed as 

^ H s^^ii tkr vi.v k ,)r < i8 » 

F 

Here k° F are points on the Fermi surface satisfying the condition q • v k o = and t k o is a 
tangential unit vector in k° F . Since the velocity v k o is perpendicular to the Fermi surface 
in k° F and also perpendicular to q, the Fermi surface has to be parallel to q in k*p. For 
convex reflection symmetric Fermi surfaces in two dimensions there are two such points 
for each given q, which are antipodal to each other. 
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Motivated by the RPA result, but envisaging already corrections beyond RPA, we 
parametrize the dynamical Stoner factor for small q and u, with v/\q\ also small, as 
follows | 

Within RPA, the length scales £o an d £ and the velocity w are related in a simple manner 
to the expansion coefficients of 11° and the coupling function g(q). Assuming that the 
latter can be expanded as g(q) = g + (felql 2 + • • • for small q, we get 

e = -gc(T)-g 2 a(T) , (20) 

(WO 2 = S^^l-gaiT), (21) 

n(q) = ±— . (22) 

For g < the velocity w(q) is always positive.— The static Stoner factor S d diverges 
at the Pomeranchuk transition, if it is continuous, and the correlation length £ diverges 
accordingly as \J~Sd~. The relation for £ makes sense only if the right hand side is positive. 
This is guaranteed if we restrict ourselves to systems where the Pomeranchuk transition is 
the leading instability. For — g c(T) — g 2 a(T) < a charge density wave instability with 
a wave vector q ^ sets in first. The parameters £o and w(q) remain finite at the Pomer- 
anchuk transition and do not vary much in its vicinity. The correlation length £ (5, T) near 
the transition depends sensitively on control parameters 5, such as the chemical potential, 
and on the temperature. If the transition is continuous, £(5, T) diverges for T — > T C (S). 
Within RPA, f(5,T) diverges as (T - T c )- 1 / 2 if T c {8) > 0, and as T" 1 if 5 is tuned to a 
quantum critical point S c . 

For small q and u, with u/\q\ also small, the d-wave density correlation function can 
be written as iV^q, v) = — K° d Sd(q, v), where n° d = — lim q ^ lim^o n^(q, v) is the static 
d-wave compressibility— Its imaginary part is then given by 

V 

ImAWq, v) = -«S , (23) 

'[(&/«)' + Slq|f + (;&)' 

For £ 3> £o an d small finite q this function exhibits a pronounced peak at low frequencies 
with a steep slope for v — ■> 0, in agreement with recent numerical results for iV d (q, v) 
obtained for the t-J model within slave boson RPA.— 

We now discuss corrections due to contributions beyond RPA. The exact density cor- 
relation function iVd(q, v) can be written in the form of Eq. (jSJ, with the full polarization 
function Ild(q, u), which is dressed by interactions, instead of the bare one. Analogously 
the full effective interaction is given by Eq. JHJ) with dressed polarization functions II (q, v) 
and n d (q, v). Close to a continuous phase transition two types of interaction corrections 
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can be distinguished, namely regular interactions which remain finite at the transition and 
singular effective interactions associated with large order parameter fluctuations which di- 
verge. 

Corrections to RPA and corresponding subleading corrections to Fermi liquid behavior 
due to regular interactions, in a generic stable Fermi liquid regime, have been analyzed 
thoroughly in the last few years.-^ The low energy behavior of most quantities receives 
non-analytic corrections to Fermi liquid behavior in dimensions d < 3. For example, in 
two dimensions the spin susceptibility varies as |q| instead of |q| 2 for small q at T = 0. 
Remarkably, the charge susceptibility for small q remains unaffected. In this case non- 
analytic contributions appearing on the level of single Feynman diagrams cancel sys- 
tematically when all relevant diagrams at a certain order are summed.^ The arguments 
establishing the cancellation of non-analytic corrections for the charge susceptibility can 
be readily extended to our case of a d-wave density instead of the conventional density 
operator .^2 The crucial point is that a perturbing field coupling to the d-wave density op- 
erator does not alter the singularities in the polarization function at q = and q = 2kp. 
Hence, corrections due to regular interactions may shift the parameters with respect to 
the RPA result for N d (q, v) and r kk /(q, v\ but they do not yield any qualitative changes. 
Within the f-model, one such correction appears already on mean-field level: the w-term 
in ./kk'(ci) generates a constant (momentum- independent) Hartree self-energy correction, 
which renormalizes the relation between [i and density- 
Near a continuous phase transition order parameter correlations are strongly renormal- 
ized with respect to the mean-field or RPA result, if the dimensionality of the system is 
below the upper critical dimension.— These renormalizations are most naturally treated 
by a renormalization group analysis of an effective field theory, where the order parameter 
fluctuations are represented by a bosonic field. The propagator of that field corresponds 
to the order parameter correlation function, that is to A^(q, v) in our case. The singular 
effective interaction between electrons, which is generated by large order parameter fluc- 
tuations, is then mediated by the bosonic field. In the bosonic representation, corrections 
to the RPA result for the order parameter correlations are due to interactions of the Bose 
field. Close to a continuous Pomeranchuk transition at finite T c these terms are relevant 
and lead to the classical non-Gaussian asymptotic behavior of the Ising universality class 
in two dimensions. In the following we focus however on the behavior in the quantum 
critical regime near the zero temperature critical point at S = S c . In that regime the 
upper critical dimension separating Gaussian from non-Gaussian behavior is d c = 4 — z, 
where z is the dynamical exponent.— In our case of a charge instability at q = one 
has z = 3, in complete analogy to the ferromagnetic quantum critical point in itinerant 
electron systems.— Hence we have d c = 1, while the dimensionality of our system is two, 
such that Gaussian behavior is stable. However, as first pointed out by Millis,— the ir- 
relevant quartic interaction of the order parameter fluctuations changes the temperature 
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dependence of the correlation length near the quantum critical point completely compared 
to the RPA result. In particular, in a two-dimensional system with z = 3, the correlation 
length at S c behaves as^ 

({S " T) -vfkf\ (24) 

instead of the naive T~ ^divergence. 

In summary, the d-wave density correlations and the singular part of the effective in- 
teraction in the quantum critical regime can be parametrized by the dynamical Stoner 
factor Sd(q, v) as in Eq. (fT^j). where the parameters £ an d w(q) do not vary much, while 
the correlation length £ (8, T) depends sensitively on control parameters 8 and tempera- 
ture. The asymptotic behavior of £(5, T) in the critical region is strongly influenced by 
interactions of order parameter fluctuations. 



IV. SELF ENERGY 

We now analyze the low energy behavior of the self-energy in the presence of strong d- 
wave Fermi surface fluctuations. We compute the self-energy to first order in the singular 
interaction Fkk'(q; v \ non-selfconsistently as well as selfconsistently, and then discuss the 
role of vertex corrections. 



A. Random phase approximation 

To first order in T, the self-energy is given by the Fock diagram in Fig. 2. The Hartree 
term vanishes because the expectation value of the d-wave density operator vanishes in 
the symmetric phase. The Fock diagram yields 

E(k, ico n ) = ~ T Y.f 702 r kk(q, G(k + q, ico n + w n ) e * 0+ (^) , (25) 

where G is the propagator of the interacting system in a self-consistent perturbation 
expansion, which is replaced by the bare propagator Go in the non-selfconsistent version. 
Note that we have approximated Tkk' with k' = k + q by Tkk , which makes almost no 
difference since only small q contribute and the effective interaction does not vary rapidly 
as a function of k and k'. Analytic continuation to the real frequency axis yields 

£(k,cu + *0 + ) = -~ J dv J [6(1/) Imr kk (q, v + *0+) G(k + q, v + to + i0 + ) 

- f(u) r kk (q, u-u- i0 + ) ImG(k + q, v + iO + )l , (26) 



where b{y) = [e f3u — 1] 1 and f{y) = [e^ u + 1] 1 are the Bose and Fermi functions, 
respectively. It is convenient to focus on the imaginary part of E, from which the real 
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part can be easily recovered via the Kramers- Kronig relation. Using Eq. (9) to express 
the effective interaction by the dynamical Stoner factor, ImE can be written as 

Im£(k,u;) = -^ JduJ -0^[b(u) + f(p + Lu)]lmS d (q,u)lmG(k + q,uj + u). (27) 

Here and in the following G, S, and are retarded functions, that is the real frequency 
axis is approached from above. The imaginary part of Sd can be written as [see Eq. (fT9"|) ] 

ImSi(q ' v) ~ *» + [(&/«)' + (&|q|)1 ! 'Mq)|qlP ' ( ' 

In the non-selfconsistent calculation one can use ImG (p, u) = —tt5(uj — £ p ), where £ p = 
e p — yU, to perform the ^-integral analytically, which yields 

f d 2 q 

ImS(k, u) =gdl J [&(f k+q - w) + /(&+ q )] Im5 rf(q> fk+q - w) . (29) 

Since we are interested in the renormalization of low energy excitations, with k close 
to the Fermi surface, and since only small momentum transfers q contribute to the self- 
energy, it is convenient to introduce a local coordinate system in momentum space, cen- 
tered around the Fermi point kp which is reached from k by a normal projection (see 
Fig. 3), such that the vector k — kp is perpendicular to the Fermi surface in kp. We 
can then parametrize k by the variable k r = ±|k — kp|, with a positive sign for k on the 
exterior side of the Fermi surface, and minus inside. The momentum transfer q can be 
parametrized by a radial variable q r and a tangential variable q t , as shown in Fig. 3. 

The excitation energy £ k+q appearing in the above expressions for ImE can be expanded 
as £k+ q = v-k F k r + Vk F q r + for small q and k near kp. The parameter m t is given 

by the second derivative of £k i n tangential direction, m~[ l = cf t £k|k=k F - The term of 
order q\ has been included since some asymptotic results are dominated by contributions 
with | (ft | 3> \q r \- It is convenient to use q' r = q r + <jf instead of q r as integration 

variable (in addition to q t ), since the excitation energy £k+q = vy_ F {k T + q' T ) is linear in 
that variable; the corresponding Jacobi determinant is one. 



1. Ground state 

At T = the combination of Bose and Fermi functions contributing to the RPA self- 
energy reduces to 

... , f - Q(-u < v < 0) for u > , or . 

b(p)+f(u + u) = \ (30) 
I 6(0 < v < —uj) for uj < 

where 6(.) = 1, if the inequalities in the argument are satisfied, and 6(.) = otherwise. 
In the following we restrict to the case u > in derivations for definiteness, but state 
final results also for to < 0. 
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At the quantum critical point (T = 0, £ = oo), the non-selfconsistent RPA self-energy 
Eq. (|29j) can be written as 



dq' r [ dq t u(q) |q| [w - f ki ,(£v + q' r )) 



/"if 
-kr 2tt 7 2tt [a; - v^ F (k r + q' r )f + £ 4 Kq)] 2 |q| 6 1 J 

for oo > 0. One may impose a cutoff on the ^-integral, which however does not affect the 
asymptotic behavior for small uj. 

For k = kp, that is k r = 0, the asymptotic (^-dependence of ImE can be extracted by 
introducing dimensionless variables q r and q t , which are defined by q' r = {jjo/v-k F )q r and 
Qt — (£o M k F )~ 1 ^ 3 ^ 1//3 %, respectively, where Uk F = u(ty F ). For small to and k r = the 
above q-integral is dominated by almost tangential q- vectors, that is \q t \ ^> \q r \; more 
precisely \q' r \ scales as \qt\ 3 , and \q r \ consequently as \qt\ 2 - Hence, we can replace |q| by 
\q t \ and u(q) by w kF . This yields, for u> — > 0, 

t w , v ^k F % 7 > 2/3 [°°dq t \qt\(l-q r ) 
ImS(k F ,^)^— -£ F / — / — . (32) 



WkF £}/3 y 2tt y_ oc 2vr (1 - g r ) 2 + $ 

This asymptotic result does not depend on any cutoff. The definite integral can be done 
analytically; the result is (4\/37r)~ 1 . For u < one obtains the same with (— uj) 2 / 3 . Hence 
we have shown that 

Q d 2 ZL^^ 

ImS(k F , uj) " kF M 2 / 3 (33) 

for small |a>|. Note that v = V\^ F q' r — uj = (q r — l)u vanishes faster than |q| for uj — > 0, 
which justifies our expansion of S^q, z/) for small ratios ^/|q|. 

Not unexpectedly, Im£(kp,u;) has the same energy dependence as for the quantum 
critical points near phase separation 23 and ferromagnetism^i in two dimensions, and also 
for fermions coupled to a [/(l)-gauge field . 2 ^ 26 In both cases the fluctuation propagator 
has the same singularity structure as our Stoner factor S^q, v) for £ = oo. Different is, 
however, the d-wave form factor making Im£(kp,u;) strongly anisotropic. It is strongest 
near the van Hove points, while the leading terms vanish on the diagonal of the Brillouin 
zone. Subleading terms and contributions from interactions with large momentum trans- 
fers generate at least conventional Fermi liquid decay rates (of order cu 2 log |u;|) on the 
diagonal, but faster decay may be obtained due to higher order processes which couple 
different parts of the Fermi surface. 

A strongly anisotropic decay rate for single-particle excitations following a power-law 
with exponent 2/3 has also been found for an isotropic continuum (not lattice) version of 
our model. 6 However, that result was obtained for the symmetry-broken "nematic" phase, 
and the large anisotropic decay rate is due to the anisotropy of the nematic state and its 
Goldstone modes. At the quantum critical point the decay rate of the isotropic system 
also obeys a power law with exponent 2/3, but now with a constant prefactor over the 
whole Fermi surface. 



12 



For k 7^ kp, that is finite k r , the asymptotic behavior of ImX(k, uj), Eq. (|31|). is obtained 
by rewriting the integral with a dimensionless variable q r defined by q' r + k r = (u/vk F ) q r , 
and q"t defined by qt = (^qw) -1 /^ 1 / 3 ^. Here we approximate w(q) by a constant w for 
simplicity In the limit uj — > one can replace |q| by \/qf + k%. Carrying out the q r - 
integral one then obtains (including the case uj < 0) 



Im£(k, w) = jj- / — y/qf + k 2 log 

27rv kF j o 2 7T 



i^ 1 



(q 2 + K 



2\3 



(34) 



ImS(k,c) = -V^-TT7^- (36) 



with k = {u^/\uj\) 1 ^ k r . Two different types of asymptotic behavior are separated by the 
frequency scale 

^k r = <o %■? ■ (35) 

For | a; | ^> u kr , corresponding to k < 1, one recovers the result Eq. (|3^j). while for 
\uj\ <C Uk r one can expand the integrand in k~ 1 and obtains 

9dj F 1 
67T 2 W kF u^kf 

In the latter limit momentum transfers normal to the Fermi surface dominate, and the 
above result can thus be easily generalized to a direction dependent w(q) replacing u by 
■u(nk F ). Note that for small k r , the low frequency behavior of ImX(k, uj) deviates from 
|c<j| 2 / 3 -behavior only below a very small scale of order \k r \ 3 . The same crossover behavior 
has already been obtained previously for fermions coupled to a gauge field.— 

We now analyze the behavior of Im£(kF, uj) at T = in the symmetric phase at some 
small distance from the quantum critical point, where the correlation length £ is large, but 
not infinite. Of interest is in particular at which scale the |a;| 2 / 3 -behavior of Im£(kp, uj) is 
affected by a large finite £. Introducing the same dimensionless variables q r and q~ t as for 
£ = oo, one finds that the q-integral is still dominated by tangential momentum transfers 
for small uj, such that we can approximate |q| by \q t \ and w(q) by «k F - The g r -integral 
can then be carried out analytically, which yields 



gdl % /3 M 2/3 f°°dq t 
Im£(k F ,w) = 7^ E ~m — / 7J-?* lo § 



1 + ^773 



2\2 



(37) 



where ( = (wk^o/M) 1 ^ -1 - ^he ^-integral can be done analytically, leading however 
to a rather lengthy expression. There are two asymptotic regimes, separated by the 
characteristic frequency scale 

^ = u kF e /e = u^o's-^ 2 . (38) 

For | a; | ^> uj^ one has ( <C 1, which leads back to the result Eq. (J3~3~j) . For \uj\ <C uj%, an 
expansion of the integral yields the asymptotic behavior 

Im£(k F , uj) = ? log ^ . (39) 
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Hence, below the scale ui^ one recovers Fermi liquid behavior. Close to the quantum critical 
point, that is for large £ and Sd, the crossover from \cu\ 2 ^ 3 scaling to Fermi liquid behavior 
sets in only for very small frequencies, and the coefficient in front of the asymptotic 
uj 2 log r| law is anomalously large. A similar non-Fermi to Fermi liquid crossover occurs 
near antiferromagnetio2£ and ferromagnetic 2 ^ quantum critical points. In Fig. 4 we show 
the non-selfconsistent RPA result for Im£(kF, uj) in the ground state for various choices 
of £, as obtained from a numerical integration of the finite £ analog of Eq. QUI J). The 
choice of parameters £o — v u. F — 1 corresponds to fixing a length and energy scale. The 
somewhat arbitrary choice of u^ F is not very critical due to the weak dependence of the 
results on Uk F . 



2. Low finite temperature 

At T > the self-energy in the non-selfconsistent RPA approximation, Eq. (|29j). can 
be written as 

ImS(k,w) = gdl p J ^ J ^ {b[v kF (k r + q' r ) - uj] + f[v kF (k r + q' r )]} 

M (q) |q| [v* F (K + <£) - U(]] 

r_. n, i ~/\ . .12 i \(t- i it i„i\2i2 i n 5» ■ \ J 



To tackle the asymptotic behavior of this integral for low T, small uj and small k r , it is 
instructive to consider first the special case k r = uj = 0, that is 

ImS(k F ,0) = 9 d l F J ^ J [K^fQt) + f( v k F q' r 

K F g;) 2 + [(WO 2 + (e |q|) 2 ] 2 Kq)|q|] 2 ' 1 ] 



We introduce dimensionless variables q r and q t via the relations q' r = £ 2 q r /£ 3 and g 4 = 
respectively. We now assume that £ diverges faster than T -1 / 3 for T — > 0, as is indeed 
the case when we approach the quantum critical point from the quantum critical region.— 
Then the above integral is dominated by momenta with a small ratio v^q^/T, such that 
the Bose function can be expanded as b(vk F q , r ) — > T/(vk F q' r ), and the Fermi function can 
be neglected. Furthermore we can exploit that the integral is dominated by momenta q 
which are almost tangential to the Fermi surface for large £, since \q' r \ scales as \qt\ 3 , and 
hence \q r \ as \qt\ 2 - We can thus replace |q| by \q t \ and w(q) by u^ F . We then obtain the 
simple result 

t y /i n x 9dj F f°° dq r f°° dq t n kf \q t \ g dj F 

K F) % * loo 2vr 2vr ( UkF g r )2 + (1 + q 2 f (u kF q t ) 2 Av^ 2 5 

(42) 
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for T — > 0. A similar contribution of order T£ has been found already earlier for almost 
antiferromagneti o 2 ^ 2 ^ and almost ferromagnetic 22 metals. Note that ImE(kF, 0) does not 
obey the same power law as a function of T as Im£(kF, uj) as a function of u; at T = 0. 



The T 2 / 3 -law proposed in Ref. [l(J, which one might expect by identifying T and uj scaling, 
does not describe the leading asymptotic behavior at low T. 

We now generalize the preceding analysis to finite ou, which shall however be sufficiently 
small that we can still use the expansion of the Bose function and neglect the Fermi 
function. To this end we set uj = v^ F x/C,, where x is a dimensionless scaling variable 
which is kept fixed in the low temperature limit. Once again we introduce dimensionless 
integration variables q r and q t , now defined by q' r — uj/vk F = £ 2 g r /£ 3 and q t = qt/£,, 
respectively. For £ — > oo we can then replace q r in |q| = \/ q% + <? 2 by uj/vk F , which yields 
|q| = g/£ with q = a/ qf + x 2 . The Bose function can again be expanded and the Fermi 
function neglected for T — ► 0, if £ diverges faster than T -1 / 3 . We then obtain 

ImS(k F ,a;) -> / T£l(x) (43) 



4^k F £o 2 



with a dimensionless scaling function 



]( T ) ^ f°° dqr [°° dqt Av kF u(q)q 

{ ) J-oo 2vr 2vr (v kF q r y + (1 + g 2 ) 2 [«(q) g] 2 1 1 

The unit vector q can be parametrized by x and qt, it does not depend on q r . Doing the 
elementary ^-integral, the velocities t>k F and tt(q) in the above expression for l(x) drop 
out completely. Carrying out the remaining ^-integral, we obtain the simple universal 
result 

l(x) = -=L= . (45) 



X 



For uj = one has 1(0) = 1 and the above special result for ImX(kp, 0) is recovered. 
For large x the scaling function decays as x^ 1 . Hence, the contribution from the Bose 
function singularity to Im£(kp,ci;) leads to a peak with a height scaling as T£(T) and 
a width of order i> ki? /£(T). The product of height and width is thus proportional to T. 
Since the contribution proportional to the Bose function to Im£(k, uj), Eq. pOj). depends 
only via the linear combination uj — Vk F k r on uj and k r , the right hand side of Eq. (|4~H|) is 
applicable also for k ^ k^, where it yields the contribution from the expanded Bose term 
for T — > 0, uj — > Vk F k r , with fixed x = (uj /v^ F — k r ) £ . 

The above asymptotic result is entirely due to "classical" fluctuations, corresponding 
to the contribution with v n = to the Matsubara frequency sum in Eq. fl25(l . The analytic 
continuation of that contribution to real frequencies reads 

S c (k, u ) = -tJj0z T kk (q, 0) G(k + q, uj) (46) 
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Note that Tki^q, 0) is real and does not depend on the parameter w(q). For G = Gq one 
can easily show that indeed 

ImS c (k,^) - p^T£l[(u/v kF - kr)$ (47) 

with l(x) from Eq. (|45j) . A similar scaling behavior of the self-energy in almost antiferro- 
magnetic metals has been derived in Ref. 1^. 

We now split the total self-energy as £ = S c + S g , where the "quantum" contribution 
is obtained by summing Matsubara frequencies v n ^ in Eq. After analytical 

continuation to real frequencies, ImS c was obtained from the Bose function singularity 
T jv. To analyze ImS 9 for real frequencies we thus subtract T/u from the Bose function, 
that is we replace b{y) by the regular function b{y) = b(u) — Tjv in Eq. (JJUJ). The 
asymptotic behavior of ImS 9 at low energy scales (low frequency and temperature) can 
be extracted by using the same dimensionless variables q\ and q r as already in the case 
T = 0, and scaling uj as T by keeping uj = uj/T fixed in the limit T — * 0, uj — > 0. 
Asymptotically one can replace |q| by \qt\ and u(q) by Uk F as for T = 0. Furthermore 
one can neglect £ -2 in the denominator of ImS^ since £ -2 scales to zero faster than |q| 2 . 
The ^-integral can then be carried out analytically, and we obtain 

ad 2 u 1/3 U 2/3 

Im£«(k F , uj) - y -^- kF ' 3 ' s{u) , (48) 

with the universal dimensionless scaling function 

1 1 1 



sgn(^) /"°° dg r 



3^3 7-oo 2vr 



; o)(gr-l) _ i ^(q? _ 1) e ^9r + i 



g,| 4 / 3 ' 



(49) 



Note that the above integral converges since the Bose function pole has been subtracted. 
A plot of s(u) is shown in Fig. 5. For |cD| — ^ oo the scaling function tends to 4 ^ , and one 
recovers the zero temperature result, Eq. (}3*3*j) . The convergence to the zero temperature 
limit is however rather slow. For small \uj\, s(u) is negative and proportional to |u>|~ 2 / 3 , 
such that 

Im£ 9 (k F ,0)^a^-^-^T 2 / 3 , (50) 
v kF £ 4 / 3 

where a ~ —0.15 is a numerical constant. Note that ImX 9 (kF, 0) is positive but smaller 
than the absolute value of the classical contribution for low T, since the latter is propor- 
tional to T£, such that the imaginary part of S c + Tfl remains negative, as it should. For 
k 7^ kp, that is for finite k r , the momentum dependence of Im£ 9 (k, uj) is negligible for 
\uj\ ojk r , with L)k r = u^\k r \ 3 , as in the zero temperature case. 

To summarize, at low finite T the RPA self-energy is given by a classical contribution 
S c of order T£, see Eq. (J47|) . and a quantum contribution S 9 of order T 2//3 and leu] 2 / 3 , 
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which obeys u/T- scaling. For |u;| ^> u>k r the latter is described by Eq. (|48p. A similar 
structure of the self-energy, with a classical part and a quantum part obeying (uj/T)- 
scaling, has been obtained already earlier for electrons coupled to strong ferromagnetic 30 
or antiferromagnetic 27 fluctuations in the quantum critical regime. In the latter case the 
self-energy is singular at special hot spots on the Fermi surface, and S 9 scales with an 
exponent 1/2 instead of 2/3. In Fig. 6 we show results for Im£(kp,u;) as obtained from 
a numerical integration of Eq. (pfUJ) at various temperatures. The correlation length has 
been chosen as £(T) oc (T| logT|) -1 / 2 , that is the temperature dependence derived by 
Millis. 22 Analo gous results for k ^ kp are plotted in Fig. 7. In agreement with the above 
asymptotic analysis the graphs show a dispersing finite T structure sitting on top of an 
almost momentum independent background which is proportional to \cu\ 2 ^ 3 for T <C \u)\- 

B. Self-consistency 

The self-energy obtained above modifies the propagator G strongly at low energy scales. 
We now check the self-consistency of the results obtained in the preceding subsection, 
that is we analyze to what extent the self-energy obtained from the self-consistent RPA, 
Eq. (|2*5|) with a dressed propagator G, differs from the one computed with G . 

1. Ground state 

At T = and for u > 0, ImS(k, to) in self-consistent RPA, Eq. (12 7)1 . can be written as 

Im£(k,u;) = -^ /" ' de f ^ ImS^q, e - u) ImG(k + q, e) . (51) 

7r J J 2n J 2n 

We focus on the quantum critical point, £ = oo, such that 



ImS^q, v) = 



u(q) 


q| 


V 




(q)] 2 |ql 


6 



(52) 



Motivated by the perturbative results in Sec. IV. A, we assume that the low-frequency be- 
havior of Im£(kF, u) is given by C2|w| 2//3 , where Gi is a negative constant, and ImS(k, uj) 
with a small distance k r from the Fermi surface is given by the same behavior for frequen- 
cies above the small scale Uk r = u^q k 3 . The analytical properties of S in the complex 
plane then dictate 

..poo ft I /■ 1 2/3 

ReS(k F ,cu) = --V \ did' 2 1 ' =C 1 sgn(cu)|cu| 2 / 3 , (53) 

where V denotes the principal value and C\ = v^CV Anticipating that for small u and 
k = kp the integral in Eq. (joTjl is dominated by contributions with e ^> u^qf, we can 
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approximate 

lmG(k F + q, e) = 2 (54) 

[e - v^ F q' r - Cisgn(e)|e| 2 /' 3 ] 2 + {C 2 \e\ 1/3 Y 

under the integral. Scaling out the u;-dependence from the integral in Eq. (J5TJ) one finds 
that the integration variable e scales as lo, qt as uj 1 / 3 , and q' r as uj 2 /' 3 . Hence, for small uj 
one can replace |q| by \q t \ and -u(q) by u^. This implies that the ^-integration acts only 
on ImG, yielding simply 

/f , mG(kF + q , f ) = -^_. (55) 

This is independent of C\ and C 2 , which means that the self-energy drops out completely! 
Carrying out the integral over e and qt one then recovers the result for ImX(kF, u) obtained 
already within the non-selfconsistent RPA. Hence, the replacement of Go by G does 
not affect the asymptotic low frequency behavior of ImS(ki?,u;) at all, it does not even 
modify the prefactor. The same result has already been obtained earlier in the formally 
similar problem of fermions coupled to a U (l)-gauge field,— and also for antiferromagnetic 
quantum critical points.— 

The above arguments can be easily extended to the case k ^ kp with a small finite k r , 
since the latter affects ImE(k, uj) only for frequencies below the small scale Uk r oc k 3 . 



2. Low finite temperature 

We first analyze to what extent the peak-shaped classical contribution to ImS(k, uj) 
is modified by self-energy feedback into G. As before, we compute this contribution by 
expanding the Bose function, b(u) « T/v, which yields 

ImS c (k, u ) = -i^T f — [^ lmS d (q, v) ImG(k + q, oj + u) , (56) 

where G is the full propagator. To scale out £ from ImS^q, z/), we use dimensionless 
variables v and q defined by v = t>k F £o ^ an d |q| — respectively, such that 

t q ( x £ 2 i; kF »(q)gi> 

ImS d (q, v) = -J 2 , 2 , ~ 2 , 2r / M -i 2 • (57) 

The inverse propagator can be written as G" 1 (k + q, uj + v) = uj + v — v kF (k r + q' r ) — 
E(k + q, uj + v). 

Consider first the case k = kp and uj = 0. Then S(k + q, uj + v) can be replaced by 
£(kp + q, 0) for T — > 0, since v scales to zero as £~ 3 . Motivated by the perturbative 
result we assume that ImS(k, 0) is of order T£ for k = kp and decreases for momenta 
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away from kp. Note that T£ oc (^log^)" 1 , since £(T) oc (TlogT) -1 / 2 . Then v can be 
neglected completely in the expression for the propagator, such that G _1 (kp + q, v) = 
—Vk F q' r — E kF (g£,, 0), where E kF (g^,, to) = E(kp + q, to). Now the integral over v acts only 
on ImSd, yielding 

[dp £ 2 1 

/ — lmS d (q,u) = 7T-2 — — 2 = 7r5d(q,0) . (58) 

Writing g 2 = g 2 + g 2 , we can also carry out the integral over g t and obtain 

dq t [ dv vr £ 1 

— / — ImS d (q, ^) = --2 y= — = • (59) 

Using q' r ~ q r = q r /£, and collecting all terms we get 

ImS c (k F , 0) = T e / — ; 1 Im ^— ^ , . (60) 

The perturbative result for the self-energy indicates that £ kF (g r /£, 0) is °f order T£. For 
£(T) oc (T\ logTl)" 1 / 2 the self-energy correction ££ kF (g r /£, 0) in the above integral is then 
suppressed as l/log£ for T — > (£ — >■ oo) and can thus be neglected compared to q r , 
albeit only with logarithmic accuracy, such that Im[-u kF g r + £E kF (g r /£, 0)] _1 — > 7ru k <5(gv) 
and we recover the perturbative result Eq. (pf2*j) for lm£(kp,0). More generally, the 
perturbative result is not modified (asymptotically) in a self-consistent treatment as long 
as £(T) increases slower than T _1,/2 | logTl 1 / 2 for T — > 0.— 

For general k (at or away from kp) and finite uj we can still neglect v in G, such that 
G _1 (k + q, oj + v) = uj — v-k F (k r + q' r ) — T l ) iF (k r + q l r ,u}). Then ImS^ can again be integrated 
independently over v and q t , which yields 



2£o V 27r ^/l + g2 t; kF g r + £[?; kF A; r - to + S kF (£v + q r ' 

(61) 

The classical contribution to the self-energy feedback in the above integral is again sup- 
pressed at least as l/log£, such that only the quantum contribution S 9 remains on the 
right hand side. The latter is not negligible at u ^ 0, it rather dominates over the bare fre- 
quency dependence of G for small finite uj. Let us assume that the RPA result for Tfl is not 
modified qualitatively by self-consistency, as will be indeed verified below. For k r = or k r 
scaling to zero more rapidly than lu] 1 ^ 3 one thus has simply £ kf ,(fc r + q r /£, uj) — > £ 



uj 



where E k (u;) is of order \u\ 2 ^ 3 for small uj,T with T <C \u\ and of order T 2//3 for small 
uj,T with \u\ <C T. In the latter case Tfl is smaller than E c , such that the self-energy 
feedback can be neglected completely. For \ui — Vk F k r \ "C |S kj? (o;)|, the term — v\^ F k r 
can be neglected in Eq. (}6T| . such that Im£ c (k, a;) becomes independent of k r and the 
frequency dependence is given by a scaling function which decreases monotonically with 
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a dimensionless scaling variable proportional to £|cj| 2 / 3 . Hence, in the self-consistent cal- 
culation Im£ c (k, £k) is not fc r -independent but rather decreases with increasing \k r \ on 
a scale of order £~ 3//2 , that is quite rapidly. For \u — v^ F k r \ ^> |££ (u;)| the self-energy 
feedback is negligible and one recovers the non-self consistent RPA result. In general one 
will find a crossover between the limiting cases. In particular, the decrease of Im£ c (kF, u) 
as a function of increasing \u\ occurs on a scale of order £ _3//2 for the smallest frequencies 
and then, more slowly, on the larger scale For k 7^ kp the self-consistent result for 
Im£ c (k, uj) depends on k r and u in a more complicated way than just via the difference 
^ — Vk F k r . Asymptotically it depends only on to if k r scales to zero faster than |a>| 2 / 3 . 

We now check possible modifications of S 9 due to self-energy feedback at low finite 
temperatures. Subtracting the pole from the Bose function, the quantum contribution to 
Eq. (J27)) can be written as 

ImE , (k , ^ ) = -iAj dv JKJ d Ji m _ I + /( „ + a)] 

xlm5 d (q, v) ImG(k + q, to + v) . (62) 

We consider the case k = kp. Frequencies are scaled with temperature, that is u = Too 
and v = TP. At low T and uo one can, once again, replace |q| by \qt\ and u(q) by Wk F 
in lmSd(<l,v), and neglect £~ 2 . Then the g^.-integration acts only on ImG. The inverse 
propagator can be written as G~ 1 (ki? + q, u + u) = uj + v — v-^ F ^ r — T,l_ F (q' r ,LJ + v) — 
T!^ F {q' r) uj + v). The question now is whether H^ F (q' r ,tj + v) leads to a significant q' r - 
dependence, which could spoil the simple result Eq. for the ^-integral of ImG. Since 
uj + v scales to zero as T, and thus faster than £~ 3 / 2 , one can replace T,^ F (q' r ,uj + v) by 
(q' r ,0). The latter tends to the constant (0, 0) for \q' r \ <^ The largest \q' r \ 
contributing to J dq' r ImG are proportional to T£, which indeed scales to zero faster than 
for T — > 0, albeit only by a factor of order logT. Hence, to logarithmic accuracy, 
we can indeed neglect the (^-dependence of the self-energy in f dq' r ImG, such that the 
simple formula Eq. (jHHj) is still valid, and the non-selfconsistent result for Im£ IJ (kF, uS) is 
confirmed. The same result is obtained by the same arguments for k 7^ kp. 

A comparison of the asymptotic result for S c with a numerical evaluation of the self- 
consistency equations reveals that the feedback of S c into Eq. (jfiT|) becomes negligible 
only on extremely small energy scales. This is not surprising as we have shown that this 
feedback vanishes only logarithmically. By contrast, the effects of self-energy feedback 
into S 9 are indeed relatively small also for finite temperatures and frequencies. In the 
following, we therefore use the non-selfconsistent RPA result for Yfl, but compute S c from 
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a self-consistent solution of the complex version of Eq. (|61jh that is 




1 



(63) 



Vk P (q r + K) - u + T^ F {k T + q T ,u) + ££ p (A; r + q r ,u) ' 



This equation can be easily solved numerically by iteration. Due to the weak radial 
momentum dependence of S 9 one can neglect the g r -dependence in its argument. Instead 
of using the asymptotic scaling form of S 9 , which is valid only at sufficiently low uo 
and T, we compute S 9 by integrating Eq. fiT)|) numerically and subtracting the classical 
contribution. In Figs. 8 and 9 we show results for Im£(k, u) as obtained by the above 
procedure. One can see that the self-energy feedback suppresses the peak generated by 
classical fluctuations, the effect being of course stronger for larger gd^ . 

In summary, self-energy feedback in a self-consistent calculation affects only the con- 
tribution from classical fluctuations to E(k, u) at T > 0. 



At zero temperature, vertex corrections and their feedback on the self-energy have been 
analyzed in detail by Altshuler et alM for fermions coupled to a £7(l)-gauge field, which 
share many features with our model. They showed that vertex corrections may lead 
to moderate finite renormalizations, but the qualitative behavior of the self-energy, in 
particular the power- law with the exponent 2/3, remains unchanged. This was confirmed 
by a renormalization group analysis of the gauge theory, and also for the formally similar 
case of a quantum critical point near phase separation.^ The arguments used in the above 
works can be directly transferred to our system and will not be repeated here. 

We have not performed a complete analysis of vertex corrections at finite temper- 
atures. By virtue of cu/T-scaling, one may expect that quantum contributions (from 
finite Matsubara frequencies) to vertex corrections behave similarly at zero and low finite 
temperatures, and lead to finite renormalizations only. By contrast, contributions from 
classical fluctuations at T > have no counterpart at T = and may thus behave dif- 
ferently. We have therefore analyzed the first order vertex correction (see Fig. 10) due to 
classical fluctuations in the quantum critical regime. In the static limit, v = and q — > 0, 
the classical contribution to the first order vertex correction is given by 



C. 



Vertex corrections 




(65) 



(64) 
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For T — > 0, the self-energy term £Ek F (g r /£,0) is of order T£ 2 for small q r , and vanishes 
thus logarithmically. The above integral diverges thus logarithmically as (T£ 2 ) -1 , which 
cancels precisely the prefactor. Hence, the classical vertex correction remains finite. 

We have not analyzed any higher order vertex corrections. However, by virtue of the 
above results and arguments we do not expect that vertex corrections modify the RPA 
results for the self-energy drastically in the quantum critical regime. 



V. SINGLE PARTICLE EXCITATIONS 



The momentum resolved spectral function for single particle excitations is given by 

1 1 _ 7T- 1 |ImS(k,cj)| 

( ,U >~~n m u-&-Z(k,tj) ~ [uj - & - Re£(k, uj)] 2 + [ImS(k, uj)] 2 ' ((>(>) 

where = e k — ^ an d E(k, uj) is the self-energy computed in the preceding section. 
Close to the Fermi surface the self-energy for our model obeys the symmetry relations 
ReSk F (— k r , — uj) = — ReSk F (fc r , ui) and Im£k F (— k r , — uj) = ImSk F (fc r ,^), which implies 
Ak F (—k r , —uj) = Ak F (k r ,ui), since £k = Vk F k r for small k r . 



A. Ground state 

In Fig. 11 we show results for the spectral function A(k, uj) as a function of uj for several 
choices of k r . The underlying self-energy has been computed at the quantum critical point 
(T = and £ = oo) by integrating Eq. (|31|) numerically. In the following we discuss the 
most important features by using the analytical results from Sec. IV. 

At the quantum critical point the asymptotic low-energy result for the self-energy can 
be summarized as 

S(k,w) -> X kF M = ~ C k F 

Strictly speaking this simple fc r -independent behavior is valid only for |a?| ^> uj^, but the 
scale u)k r is proportional to k% and thus tiny for k near the Fermi surface. The prefactor 
C*k F depends strongly on the position of kp. It decreases rapidly for k^ near the Brillouin 
zone diagonal, where c?k F vanishes, while it is enhanced near the van Hove points, where 
fk F becomes small. The kp-dependence of u\[^ is comparatively weak. The competition 
between uj and the self-energy in the denominator of A(k,uj), Eq. (|66p. leads to the 
characteristic energy scale 

= Cl F oc % , (68) 
which obviously varies very strongly with kp. 



sgn(cj) + 



Vs. 



2/3 



where C\ 



dt 

k F 



U 



47TV kf 



.4/3 
Cn 



(67) 
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For fixed k with |£k| ^> uj c the spectral function A(k,u) has almost Lorentzian shape 
as a function of uj, with a maximum near uj = £k and a width of order Ck F |£k| 2 ^ 3 - The life- 
time broadening thus decreases more slowly than the excitation energy £k, as k approaches 
the Fermi surface, such that no well-defined quasi-particles exist. For |£k| « uj c the 
maximum of A(k, is shifted strongly away from £k and the width is of order of the 
peak energy. For |£k| ^ one can neglect to compared to Re£(k, to) in Eq. (jHEJ), and 
A(k.,to) is now peaked at uj = £k, with the renormalized energy scale 

4 = sgn(a)(C ki !Kk|) 3/2 ocfc r 3 / 2 . (69) 

Extracting a dispersion relation from the momentum dependence of the peak position in 
A(k, to) one thus obtains a flat band with a vanishing slope near the Fermi surface. The 
width of the peak centered around £k is of order Ck F |£k| 2 ^ 3 = |£k| and thus linear in k r . 
For k = kp the spectral function diverges as |co>|~ 2//3 for uj — > 0. 

Momentum scans of A(k, uj) perpendicular to the Fermi surface at fixed uj lead to 
Lorentzian peaks centered around k r = ^— [uj + Ck F sgn(c<j)|u;| 2 / 3 ] . Some such scans are 
shown in Fig. 12 for various choices of uj. The integral 

/■|^4(k, W ) = -J— (70) 

does not depend on the self-energy. The k-integrated density of states near the Fermi level 
remains therefore unrenormalized, that is, it is determined by the bare Fermi velocity. 



B. Low finite temperature 

At low finite T the results for the self-energy can be summarized as follows. The 
total self-energy is a sum of two distinct contributions, S = S c + S 9 , where S c is due 
to classical, and Tfl due to quantum fluctuations. The quantum contribution Tfl can be 
computed from non-selfconsistent RPA and obeys (a; /T) -scaling at very low uj and T, see 
Eq. (148)) . Its dependence on k r is very weak. The classical contribution S c is affected 
more strongly by self-energy feedback, especially by feedback of Tfl. We compute S c by 
solving Eq. (jo^j) self-consistently. The classical part violates {uj/T) -scaling and depends 
significantly on k r . 

The most significant temperature effect is that ImS(k^, 0) increases quickly from zero 
to sizable finite values upon increasing T. For small T near the quantum critical point 
we obtained 

ImE(k F , 0) - /% T£ cx ^ J -r^— (71) 

both in the non-selfconsistent and self-consistent calculation. This cuts off the divergence 
of A(k.F, uj) for uj — > occurring at zero temperature and replaces it by a maximum of order 
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a/| \ogT\/T. Note that Im£(kF, 0) vanishes much slower with T than in conventional 
Fermi liquids, where one has T 2 (in 3D) or T 2 | logT| (in 2D) behavior. 

In Fig. 13 we show results for A(k, uj) as a function of u for various choices of k r at a 
fixed low temperature. The underlying self-energy has been computed by the procedure 
described at the end of Sec. IV. B. The most striking differences compared to the ground 
state results (Fig. 11) are of course seen for k near kp, where the peaks are now much 
broader. Note also the steep shoulder near uj = for the spectra with k near kp at strong 
coupling (lower panel). The complementary view, A(k, u) as a function of k r for various 
fixed uj, is shown in Fig. 14. Here the line shape resembles a Lorentzian function with a 
relatively large width. 

VI. CONCLUSION 

In summary, we have presented a detailed analysis of quantum critical fluctuations 
and their effect on single-particle excitations in a two-dimensional electron system close 
to a d-wave Pomeranchuk instability. The fluctuations can be viewed as long-wavelength 
density fluctuations with a d-wave form factor and also as d-wave shaped fluctuations 
of the Fermi surface. They lead to a strong singularity in the dynamical d-wave density 
correlation function at small momenta and frequencies, and to singular forward scattering. 
The momentum and energy dependence of the singularity is captured essentially correctly 
by a Gaussian theory (RPA); interactions of fluctuations modify only the temperature 
dependence of the correlation length £(T). 

Single-particle excitations are strongly affected by the fluctuations. We have analyzed 
the electron self-energy S(k, u) within plain and self-consistent RPA, focussing especially 
on the low-energy behavior in the quantum critical regime. The dominant contributions 
due to singular forward scattering are proportional to c? 2 ., where is a form factor with 
d-wave symmetry, such as cosk x — cosk y . For k near the Fermi surface this leads to a 
strong tangential momentum dependence of S(k, uj). The singular contributions vanish 
on the diagonal of the Brillouin zone, and have the largest amplitude near the van Hove 
points. By constrast, the momentum dependence of S(k, uj) perpendicular to the Fermi 
surface is much weaker at low temperatures. Hence, momentum scans of the spectral 
function perpendicular to the Fermi surface have almost Lorentzian line shape. 

At the quantum critical point, the real and imaginary parts of the self-energy scale as 
\tu\ 2 / 3 with energy. This leads to a complete destruction of quasi-particles near the Fermi 
surface except on the Brillouin zone diagonal, due to the prefactor d\. The dispersion 
of the maxima of the spectral function A(k, u) flattens strongly for momenta k near the 
Fermi surface away from the zone diagonal. On the other hand, the momentum integrated 
density of states is not renormalized significantly by the Fermi surface fluctuations. 

In the quantum critical regime at T > the self-energy consists of a "classical" and a 
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"quantum" part with very different dependences on T and u. The classical part, which is 
due to classical fluctuations, dominates at u = 0, where it yields a contribution propor- 
tional to T£(T) to ImS(kir,o;). The quantum part is generated by quantum fluctuations 
and obeys (a>/T)-scaling in the quantum critical regime. 

We finally discuss whether soft Fermi surfaces and critical Fermi surface fluctuations 
could play a role in cuprate superconductors. Due to the coupling of electron and lattice 
degrees of freedom a symmetry-breaking Fermi surface deformation is generally accompa- 
nied by a lattice distortion, and vice versa. Structural transitions which reduce the lattice 
symmetry of the cuprate-planes are quite frequent in cuprates. Close to a Pomeranchuk 
instability of the electronic system, electronic properties can be expected to react unusu- 
ally strongly to slight lattice distortions which break the symmetry of the electronic system 
explicitly. Such "overreact ions" of electronic properties have indeed been observed early 
on in several cuprate compounds . 33 ' 34 In particular, a slight orthorhombicity of the lattice 
structure would lead to a relatively strong orthorhombic distortion of the Fermi surface. 
Yamase and Kohno 3 ^ invoked this idea already a few years ago to explain peculiarities of 
magnetic excitations in cuprates. Recent experiments on YBCO have established a re- 
markably strong in-plane anisotropy of electronic and magnetic properties , 2£i2L2S although 
the structural anisotropy of the CuCvplanes, which is induced indirectly by the CuO- 
chains between the planes in that material, is relatively modest.— Fermi surface softening 
with d-wave symmetry due to forward scattering interactions can naturally amplify the 
effect of a weak or moderate structural anisotropy. 

Since the Pomeranchuk instability breaks the orientational symmetry of the lattice, it 
is natural to compare with the mechanisms and consequences of stripe formation, which 
has been extensively discussed in the context of cuprate superconductors.— Static stripes 
also break the translation invariance in addition to the orientational symmetry, and their 
formation requires interactions with large momentum transfers, such as antiferromagnetic 
interactions. Many experimental observations, which have been attributed to static or 
fluctuating stripes,— actually provide evidence only for a tendency to orientational, not 
translational, symmetry-breaking, and could therefore be described equally well by a 
(incipient) Pomeranchuk instability. 

Strong Fermi surface fluctuations could be at least partially responsible for the non- 
Fermi liquid behavior observed in the "strange metal" regime of cuprate superconductors 
near optimal doping. In our model calculation we have obtained a strongly anisotropic 
anomalously large decay rate for single-particle excitations and a flattening of the dis- 
persion relation near the Fermi surface away from the nodal direction. The properties of 
single-particle excitations in various cuprate compounds have been investigated in con- 
siderable detail by numerous angular resolved photoemission experiments.— Extended 
flat bands in the van Hove region have been observed by various groups already in the 
early 1990s . 43 ' 44 i 45 Large anisotropic decay rates have been extracted from the linewidth 
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of low-energy peaks in the photoemission spectra observed in optimally doped cuprates, 
using in particular momentum scans perpendicular to the Fermi surface at various fixed 
energies . 4 ^ 4 ? The line shape of these scans is almost Lorentzian, which is consistent with 
our results. However, the frequency and temperature dependence of the self-energy ex- 
tracted from the experimental raw-data differs from what we obtained from d-wave Fermi 
surface fluctuations. On the other hand, the functional form of the self-energy chosen by 
the experimentalists may not be the only way to achieve consistency with the photoemis- 
sion data. In spite of the impressive progress in this experimental technique, the accuracy 
of such subtle properties as the energy and temperature dependence of the electron self- 
energy is still limited by resolution and background problems. 

Concerning transport, an anisotropic scattering rate with nodes on the Brillouin zone 
diagonal can very naturally account for the pronounced anisotropy between the intra- 
and inter-plane mobility of charge carriers, as pointed out by Ioffe and Millis^ in their 
phenomenological " cold spot" scenario. According to their idea, the intra-plane transport 
is dominated by quasi-particles with a long life-time near the diagonal of the Brillouin 
zone, while these carriers are not available for inter-plane transport, since transverse 
hopping amplitudes vanish on the diagonal. 

A spin dependent Pomeranchuk instability was recently invoked to explain a new phase 
observed in ultrapure crystals of the layered ruthenate metal Sr 3 Ru207, 49 and also to 
account for a puzzling phase transition in URu2Si2-— For a broader comparison with 
experimental data for cuprate superconductors and other layered materials, which might 
undergo or be close to a symmetry-breaking Fermi surface deformation, it will be useful 
to compute experimentally accessible response functions in the presence of strong Fermi 
surface fluctuations. 
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APPENDIX A: BARE POLARIZATION FUNCTION 

In this appendix we derive the asymptotic expressions for the bare polarization function 
II°(q, v) defined in Eq. for small q and v. The derivation is valid for arbitrary form 
factors dp with any symmetry (not only d-wave), in particular also for the special case 
d p = 1, for which 11° reduces to the conventional polarization function II . 



26 



At zero frequency, 11° can be written as 

nS(q,o)=/ ? ^ /(w) - /(w) 4- (AD 

with £ p = e p — /i. The numerator and the denominator of the fraction under the above in- 
tegral are odd functions of q, and the integrand is thus an even function. The denominator 
can be expanded as 

1 d 3 e 

ep+q/2 - w 2 = q • v p + - £ d * d wm* + ^(Iql 5 ) (A2) 

31,32,33 J1 JI Ji 

with v p = V p e p and ji, J2, J3 each running over the two possible directions x and y (in 
two dimensions). The numerator is expanded as 

/( W) - /(W.) = q • v P /(£p) + ^ fS^iL. 'Oi'u^ + °(\a\ 5 ) • ( A3 ) 



24 CZfi, dVh d Ph d P. 

Using V p /(£ p ) = /'(£ p )v p and 



9 2 6p 



+/'«-) E s^mm,^ (A4) 

31,32,33 J J J 



one obtains 



/(w0 = / / (^ P )+^r(e P ) E 4rfo) ( q -v P ) 2 +o(| q i 4 ) 

e P+q /2 - e P -q/2 8 ££dp h dp ja 24 

(A5) 

Note that terms involving third order derivates of e p have cancelled. Inserting this expan- 
sion in Eq. (jAl|) and using the point group symmetries of the square lattice one obtains 



nS(q,o) 



d 2 p 



/'(&) + 



16 



r(e P )Ae p 



q| 
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p;^p 



q| 



rf p + 0(|q| 4 ), (A6) 



which establishes Eq. (j!2j) and the formulae for the expansion coefficients a(T) and c(T) 
presented in Sec. III. In the special case dk — 1 the term involving /" in the above equation 
can be rewritten in the form of the term proportional to /'" by a partial integration, 
yielding the simplified formula for the conventional polarization function 



n°(q,o) 



d 2 p 



q| 



+ 0(\q\ 



(A7) 



27 



We now derive the frequency and momentum dependence of 11° (q, v) for small q and v 
to leading order in q and v. In the limit q — > 0, v — > the polarization function depends 
only via the ratio s = v/\q\ and the unit vector q = q/|q| on q and v: 

^--lw? f ^ .+£\-& - (A8) 

In the low temperature limit only momenta on the Fermi surface contribute to the above 
integral, since then /'(£ P ) — > — 8{£,p)- The real part is an even function of s which tends to 
a(T) in the limit s —>■ 0, with corrections of order s 2 . The imaginary part has the simple 
form 

ImII°(q, v) - J 5? /'(£ p ) 4 s (a - v p • q) . (A9) 
For q, v — > and small s this simplifies further to ImIL5(q, v) — > — p(q, T) s with 

P(4 T) = - | ^ /'(C P ) 4 *(v p • q) . (AlO) 

At T = the integrand in Eq. (jA9|) contains two 5-functions. The two-dimensional 
momentum integral can then be carried out, yielding 

Imn°(q, V) -»■ -— V — ! =^T^ rr , (All) 

dl4 ' ^ 4tt^ k ^ k o |t k o -V k o(q-v k o)| 1 J 



where t^o is a unit vector tangential to the Fermi surface in h° F , and the sum runs over 
momenta on the Fermi surface which satisfy the equation v k o^ • q = s. The formula (|18|) 
for p(q, T) at T = follows directly. 
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FIG. 3: Decomposition of momentum transfers q in normal and tangential components relative 
to the Fermi surface in kp. 
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FIG. 4: Non-selfconsistent RPA result for the imaginary part of the self-energy Im£(kF, u), 
divided by gd^ p , as a function of u in the ground state (T = 0) for various choices of the 
correlation length £; the other relevant parameters are £o = ^k F = Uk F = 1- 
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FIG. 5: Scaling function s(x) describing the asymptotic behavior of the quantum contribution 
to ImS(ki?,u;) at low u> and T; the horizontal line indicates the asymptotic value of s(x) for 
large x. 
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FIG. 6: Non-selfconsistent RPA result for the imaginary part of the self-energy ImE(kp, u), 
divided by gd^ , as a function of u for several temperatures T with a correlation length £(T) = 
3/lTlogTj 1 / 2 ; the other relevant parameters are £o = ^k F = u = 1. 
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FIG. 8: Selfconsistent RPA results for the imaginary part of the self-energy Im£(kF,u;), di- 
vided by , as a function of u) for several temperatures T, with a correlation length 
f(T) = 3/lTlogTl 1 / 2 and £ = ^k F = u = 1- Upper panel: |#|d£ F = 1, lower panel: |#|c^ f = 4. 
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FIG. 9: Self-consistent RPA result for ImS(k, u) as in Fig. 8, but now for k r = 0.1. 
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FIG. 10: Feynman diagram representing the first order vertex correction. 
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FIG. 11: Spectral function A(k, uj) at the quantum critical point as a function of lo for k r = 
— 0.0405 n with n = 0,1,2,..., 6. Fixed parameters are £o = v k F = u = 1. Upper panel: 
\9\dl F = 1, lower P anel: \a\ d l F = 4 - 
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FIG. 12: Momentum scans of A(k,u) at the quantum critical point for u> = —0.0015 — 0.0405 n 
with n = 0, 1, 2, . . . , 6. The parameters are the same as in Fig. 11. 




FIG. 13: Spectral function A(k,u) at T = 0.003 as a function of lo for k r = -0.0405 n with 
n = 0, 1, 2, . . . , 6. Parameters and £(T) as in Fig. 8. 
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FIG. 14: Momentum scans of A(k, w) at T = 0.003 for uj = -0.0405 n with n = 0, 1, 2, . . . , 6. 
Parameters and £(T) as in Fig. 8. 



